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Abstract. Based on an analysis of a specific electron trajectory in counter- 
propagating beams, Bell & Kirk (PRL 101, 200403 (2008)) recently suggested that 
laboratory lasers may shortly be able to produce significant numbers of electron- 
positron pairs. We confirm their results using an improved treatment of nonlinear 
Compton scattering in the laser beams. Implementing an algorithm that integrates 
classical electron trajectories, we then examine a wide range of laser pulse shapes and 
polarizations. We find that counter-propagating, linearly polarized beams, with either 
aligned or crossed orientation, are likely to initiate a pair avalanche at intensities of 
approximately 10 24 Wcm -2 per beam. The same result is found by modelling one of 
the beams as a wave reflected at the surface of an overdense solid. 
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1. Introduction 

Within the next few years, powerful lasers may be able to realise intensities of 10 23 
to 10 24 Wcm -2 in the laboratory, enabling novel physical processes to be investigated 
[H, [21 E]. One of these is the prolific production of electron-positron pairs predicted to 
occur when an electron is accelerated in counter-propagating laser beams [I]. 

Pair production has already been realised in the lab. using laser intensities ~ 
10 20 W cm -2 . The technique, proposed in [5] , involves using the laser beams to accelerate 
electrons to MeV energies, and then allowing them to produce an electromagnetic 
cascade in a foil made of high-Z target material. Depending on the experimental set- 
up, positrons are created either by the trident process, in which an energetic electron 
interacts with the electrostatic field of a target nucleus, or via an intermediate gamma- 
ray (produced by bremsstrahlung) that subsequently pair-creates in the field of a nucleus 
by the Bethe-Heitler process [HI [TJ, [HI [9] . The technique operates at relatively modest 
laser intensity, but converts only a small fraction of the laser pulse energy into pairs. 

An alternative mechanism that has been intensively studied but not yet observed 
is the spontaneous creation of pairs out of the vacuum by laser beams. Schwinger [TO] 
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predicted spontaneous pair creation in an electric field that approaches the critical value 
^crit = 1-3 x 10 18 Vm" 1 , which is achieved at a laser intensity of roughly 10 29 Wcm -2 . 
Conservation of energy and momentum forbids this process in a plane-wave beam, but it 
can operate in counter-propagating beams [HI [12] and may carry interesting information 
on the beam sub-structure [13J . The effect does not have a sharp threshold, but the rate 
drops rapidly as the amplitude E of the electric field decreases, since it contains a factor 
exp (— E cAt /E). Detailed calculations predict observable consequences at intensities as 
low as 10 26 W cm -2 [TT] , but these are unlikely to be achieved within the next few years. 

A third mechanism of pair production using an intense laser beam was realised 
in an experiment at SLAC |14| . A beam of 46.6 GeV electrons was fired into a laser 
beam, where electron positron pairs were created. This was interpreted as being due to 
a two-step process: First a GeV photon was created by nonlinear Compton scattering 
of multiple laser photons by a relativistic electron. Then this photon interacted with 
multiple laser photons to create a pair. This technique is basically the same as the 
first method described above, except that the electrostatic field of the target nucleus 
is replaced by the electromagnetic field of the laser. The possibility that pairs might 
be directly produced by the trident process in the laser fields was not discussed. The 
positron yield observed in this experiment was very small. 

The mechanism for prolific pair production suggested by [4] is closely related to 
that realised at SLAC. However, instead of using a particle accelerator, the laser beams 
themselves accelerate the relativistic electrons. The entire cycle of electron acceleration 
and pair creation (either by the electromagnetic analogue of the trident process or via 
an intermediate real photon) takes place when laser beams counter-propagate in an 
under-dense plasma. At intensities ~ 10 24 Wcm -2 the number of pairs created per 
plasma electron was estimated by jl] to be about unity. Since each new electron and 
positron is, in its turn, accelerated in the laser beams, one expects an avalanche of pairs 
that should absorb a significant fraction of the laser pulse energy. This process would, 
therefore, swamp the pure vacuum effect. 

The prediction of [I] was based on the analysis of a particularly simple electron 
orbit at the magnetic node in the field of two counter-propagating, circularly polarized 
pulses, and used simplified descriptions of the physical processes. In particular, the 
photon emissivity was treated using a monochromatic approximation. In this paper we 
examine this mechanism more closely. The analysis is extended in two ways: we use 
improved approximations to the physical processes, and embed these in a numerical 
scheme that follows the trajectory of an electron in more realistic models of the laser 
fields. Vacuum fields are used to describe the laser beams, which is appropriate if 
they propagate in an underdense plasma, and we examine the trajectories of initially 
nonrelativistic electrons picked up by the beams at different points in the plasma. 

In section [2] we give details of the relevant physical processes. These all involve 
multi-photon interactions with the two laser beams. However, we show that, in the 
parameter range of interest, they can also be viewed as interactions of an electron or high- 
energy photon with uniform, static electromagnetic fields. The processes are (a) photon 
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emission by a relativistic particle. In various special cases this is known as synchrotron 
radiation, magneto-bremsstrahlung, curvature radiation, non-linear Thomson scattering 
or non-linear Compton scattering; here we use the term "synchrotron radiation" to 
mean the generic process. We will assume the static field is incorporated exactly into 
the electron propagator, in which case the process is of first order in the fine-structure 
constant a f . (b) The second-order (in cif) process of pair production by a charged particle 
in a static field. (This involves an intermediate virtual photon and is the electromagnetic 
analogue of the trident process. In the following we simply refer to it as the "trident 
process" .) (c) The first-order process of pair production by a real photon (a synchrotron 
photon, for example) in the same static field. 

In section [3] we describe our model of the laser fields, and give the classical equations 
of motion of the electron including radiation reaction. The way in which pair creation 
by the processes described in section [2] is implemented is discussed in section 13.31 Our 
results fall naturally into two parts. In section 14.11 we present calculations in counter- 
propagating, circularly polarized pulses with the same sense of rotation of the fields. 
The motivation here is to compare our results with [I] and to understand the influence 
of the improved treatment of the photon emissivity and of different choices of initial 
conditions of the electron trajectory. In section 14.21 we then examine pair production 
in more realistic pulses, with both circular and linear polarization. These pulses have 
a finite duration and the fields are assumed to be contained within a cylindrical region 
whose axis lies along the propagation direction. Our conclusions, which confirm and 
extend the predictions made in [4], are summarized in section [5j 



2. Physical processes 

2.1. The quasi- stationarity and weak- field approximations 

For an electron in a monochromatic plane wave, the photon emissivity and pair-creation 
rate is a function of only two dimensionless, Lorentz invariant parameters [T5] . These 
are the strength parameter of the wave 

a = (1) 

mCCc;i aS e r 

where E Q is the amplitude of the electric field and cji ascr its angular frequency, and the 
parameter rj that determines the importance of strong-field quantum effects: 
eh 

V — jpj W\ ( 2 ) 

where p^ is the four-momentum of the electron, F^ u is the electromagnetic field tensor, 
and |. . .| denotes the length of the four- vector. In terms of the field components Ei, B{ 
in Cartesian coordinates, with % = 1, ... 3: 

F i0 = Ei = -e ijk B k = -F"* 

and Eijk is the antisymmetric Levi-Civita symbol. 
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The interaction of a particle with an arbitrary external electromagnetic field is 
more complex. However, two approximations simplify it considerably. The first is 
that of quasi-stationarity: "instantaneous" values of the transition rates are computed 
assuming the external field is constant in time, and these values are subsequently time- 
averaged. This requires that the variation timescale of the field is long compared to the 
coherence time associated with the interaction. In the monochromatic, plane-wave case, 
the coherence time is t co h ~ (E CTit /E ) (h/mc 2 ) [15] . Therefore, the variation timescale 
of the field is long compared to the coherence time of the interaction provided 

a > 1 (3) 

If this (Lorentz invariant) condition is satisfied, the instantaneous transition 
probabilities for an electron of given p^ can be computed in the limit a — ► oo, i.e., 
for uniform, static fields. We assume that this argument can be generalized to the 
case of counter-propagating waves each of strength parameter a. A potential problem 
arises for those configurations where there exist points in space that are simultaneously 
nodes of both the electric and magnetic fields, since there the coherence length of the 
interaction becomes large. However, the processes that are important for pair production 
are confined to the regions of strong field, where the coherence length is short. 

The situation can readily be visualized in the classical picture: the coherence length 
associated with the synchrotron radiation of a relativistic electron of Lorentz factor 7, 
is the path length over which the electron is deflected by an angle I/7, which gives the 
same result as in the quantum case: t co h ~ Trie/ (|e| Eq). Variations of the accelerating 
fields on lengthscales shorter than the deflection length occur only close to simultaneous 
nodes of the E and B fields. However, only relatively low frequency radiation is emitted 
in these regions. 

The second approximation is that of weak fields. In an arbitrary, constant field, the 
transition probabilities may depend not only on 77, but also on the two Lorentz invariant 
parameters associated with the field: / = \E 2 — B 2 \ / E 2 rit and g = \E ■ B\ / E 2 rit . 
Provided, however, that 

/ < 1 and g < 1 (4) 

and 

7? 2 »Max(/,s) (5) 

one may neglect this dependence and evaluate the transition probabilities in any 
convenient field configuration that has the same value of the parameter rj [161 US]- 
Analogously, pair production by a single high-energy photon in a laser field with a 3> 1 
can be approximated as pair production in a uniform, static field, and the transition 
probability for interaction with the virtual photons of this field depends only on the 
parameter 

X=^-K\ (6) 
where = (Hu, Hk) is the photon four-momentum. 



Pair production in counter-propagating laser beams 5 

The situations we consider in this paper fulfil both the quasi-stationary requirement 
([3]) and the weak-field conditions (TjJ and ([5]). At a laser wavelength of X^ m fim, the 
strength parameter in a single, linearly polarized beam of intensity /24 x 10 24 Wcm~ 2 
is a = 855/ 2 4 A^m, so that the stationarity condition ([3]), is well satisfied in the range 
of interest (A^m = 1, I24 = 0.1 .. . 1). It is possible to configure laser beams in such a 
way as to reduce the variation timescale. For example, in a beam that is reflected from 
an overdense solid, high harmonics are generated [TTllIH]- However, unless substantial 
power is present in harmonics numbers n > 100, the stationarity approximation remains 
valid. The maximum field strength in a monochromatic, linearly polarized beam is 
E/E CT it = 2.1 x lO -3 /^/ 2 , and the field invariants / and g in (j4j) vanish in such a wave. 
In counter-propagating beams (each of intensity J24 x 10 24 Wcm -2 ), they do not vanish, 
but, although the field reaches twice the amplitude of the individual beams, one still 
has at all space-time points f,g< 1.6 x 10 _5 /24, so that PJ is satisfied. Pair production 

— 1/2 

becomes important for 77 > 0.1, which requires relativistic electrons with 7 > 50 / 2 4 • 
In this range, the weak field condition ([5]) is easily satisfied. Consequently, the transition 
probabilities for photon production and pair production by an electron (via a virtual 
photon) and by a real photon may be taken from computations performed for static, 
homogeneous magnetic fields that have the same value of the parameter 77 or x, provided 
that this static configuration also fulfils the weak-field conditions. 

In particular, choosing a situation in which an electron of Lorentz factor 7, or a 
photon of energy emc 2 , propagates normal to a constant, uniform magnetic field B, the 
relevant parameters are 77 = ^BjB cxlt and x — ( e /2) B/B crit . Since / = (B/B crit ) 2 , the 
weak-field condition (jSJ) requires us to compute these transition probabilities for a field 
B <C -Bcrit, i-e., in the limit 7, e ^> 1. This case corresponds to the classical limit for the 
electron trajectory, where the particle energy is large compared to heB/mc, and to the 
case well above the threshold at e = 2 for pair creation, in the case of the photon. These 
probabilities were calculated in the 1950 's; they are conveniently reviewed by Erber, 
[T9] , whose results we summarize in the following. Their application to arbitrary fields, 
as outlined above, and their generalization to particles of arbitrary spin is discussed by 
Baier & Katkov [16]. 



2.2. Synchrotron radiation 

The rate of production of photons by an electron of energy 7mc 2 moving normal to a 
magnetic field of strength B = bB crit is defined as 

= Vs^b^hll (7) 

dxdt h x 

where Xi defined in ([H]), describes the energy of the emitted photon. In this field 
configuration, 77 = 76, and we can multiply this equation by 7 to display the Lorentz 
invariance of each side: 

d 2 iV r-mc 2 F{n, x) 
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V3— — a f 77 (8) 



dxdt h x 




Figure 1. The quantum synchrotron function F(r),x) as a function of the 
dimensionless photon frequency \ ( see ©)• The parameter r/, defined in 
determines the importance of strong field QED effects; the energy of the radiated 
photon would equal that of the electron at logx = log(?y/2). The corresponding 
classical functions (see ([9])) are shown as dashed lines. 



For 7 ^> 1 and & < 1, (i.e., under the weak-field conditions (00) and ([5])) the function 
F(r),x) is given in equation (2.5a) of [19], and is reproduced in our notation in 



Appendix A In the classical limit, it becomes a function of 77 and \ i n t ne combination 



xh ■ 



F(v, x) -> /synch [4x/ (3*7 



as h — > 



(9) 



where / syn ch is the familiar expression for the synchrotron emissivity (summed over 
polarizations) in the Airy integral approximation 



roc 

/synch (y) = y dtK 5/3 (t) 

Jv 



The monochromatic approximation used in [1] is equivalent to the replacement 



F{r],x) -> /mono [4x/ (3r/ 2 ) 



with 



>(y) 



8tt 
V3 



* (y - yo) 



(10) 



(ii) 



(12) 



Pair production in counter-propagating laser beams 



7 



together with the choice y = 0.29 (which is where the function f sync h(y) has its 
maximum). The coefficients in front of the 5- function in (j!2p ensure that the 
power radiated by an electron is the same in the classical and the monochromatic 
appr oximat ions . 

In our problem, the most important difference between the quantum and classical 
emissivities is that the quantum emissivity is severely depleted for hard photons, when 
X ~ T) ~ 1, and vanishes for photon energies larger than that of the incoming electron: 

F(r / ,x) = 0for X >r//2 (13) 

This has an important influence on pair production, since these photons subsequently 
have the highest probability of conversion. The effect is illustrated in figure [U It is 
associated with a reduction in the (Lorentz invariant) total power P radiated by the 
electron (integrated over all emitted photons): 

P = mc 2 ^-g(r]) (14) 

3 n 

with 

The classical result is g(rj) — 1. For rj ^ 1 the lowest order quantum correction gives 

m 

g(v) « i - ^j^-v ( 16 ) 



The full expression is given in |Appendix A 



2. 3. Trident pair production via virtual photons 

Trident pair production via virtual photons is a second order process. The rate, 
computed using a Weizsacker- Williams approximation, is presented by Erber [19] . This 
approximation involves treating the electromagnetic field of the incident electron as 
a superposition of virtual photons. The computation then essentially evaluates the 
dispersion relation for these photons in the magnetized vacuum. The imaginary part 
of the refractive index gives the absorption rate, i.e., the rate of pair production. This 
procedure should provide a reasonable approximation when the weak-field conditions 
are satisfied. 

The rate of pair production for an electron of energy 'jmc 2 moving normal to a 
constant magnetic field B = bB crit is ([19], equation (4.4)): 

7 ^r = 0M ^ 2 ^(v) (it) 



The function £l(rj) is given in Appendix B For small rj it can be approximated by 



7T 5 / 2 

^)^^(3^) 1/4 exp [-8/(3?7) 1/2 
For rj 3> 1, Erber gives the asymptotic value Cl(i]) — > {tt 2 /2) In 77. 
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2-4- Pair creation by synchrotron photons 

The photon absorption coefficient in the ultra-relativistic limit required by the weak- 
field approximation is expressed by Erber p3] as an absorption probability per unit path 
length. Since we will be interested in the propagation of photons in electromagnetic fields 
that vary in space and time, we write the instantaneous absorption probability in terms 
of the differential optical depth dr traversed by the photon in an interval of time dt in 
the lab. frame: 

% = Q 'ir^ T± M < 19) 

The function T±(x) is given approximately by 

K 2 (^) 

T±(x)~0.16 1/3V3x; (20) 
X 

For small x it is proportional to exp [2/ (—3%)] and so is exponentially small, despite 
that fact that, in the weak- field approximation, the photon is always well above the 
kinematic threshold. The function peaks at x ~ 8 an d falls off to higher x as X 1 ^ 3 - 



3. Trajectories 

3.1. Laser fields 

Bell & Kirk [I] considered a particle trajectory in counter-propagating monochromatic 
waves, with circular polarization, such that the field vectors rotate together at each 
point in space. This special case is considered in detail in Section 14.11 However, such 
an approach is limited in three respects: 

(i) Real laser pulses are of finite duration, the electrons they accelerate begin and end 
their trajectories outside of the pulse trains. 

(ii) Circular polarization is a special case, that may not be the most favourable for pair 
production. 

(iii) Monochromatic waves are an idealization. In particular, if one of the counter- 
propagating waves is produced by reflection from a solid surface, it will contain 
many high-order harmonics [17] . 

We lift these limitations by considering a range of laser waveforms as follows: 

(i) To model pulses of finite duration, we multiply the monochromatic wave by an 
envelope function /(</») that contains two parameters: the duration or length L of 
the pulses (in phase units) and the thickness A (also in phase units) of the pulse 
edges: 

UW = j 

Here = 0± = z =F t is the phase of the wave, and the upper (lower) sign refers to 
the rightwards (leftwards) propagating pulse. 



1 =F tanh 



1 ± tanh 



A 



(21) 
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Figure 2. A snapshot of the electric field (normalized to unity at pulse centre) used to 
model laser pulses of finite duration according to (f!2Tj) . with L = IQtt and A = 107r/3 
(upper panel). At t = the leading edges of the pulses meet, the snapshot is taken at 
t = —50, corresponding to 26.5 fsec before interaction for a laser wavelength of 1 fim. 
The lower panel shows the model waveform used for a reflected pulse, given by (1221) . 

With flmai = 7. 



(ii) In order to avoid the stable E = nodes of the circularly polarized waves, we 
consider counter-propagating, linearly polarized waves with aligned polarization and 
with crossed polarization. In aligned polarization the electric fields of the counter- 
propagating beams are parallel, and in crossed polarization they are orthogonal. 

(iii) The easiest way to produce counter-propagating waves may be to reflect a wave 
from a solid target. To account for the complex harmonic structure of the reflected 
wave pTTJ, [18] we consider a leftward propagating wave with electric field: 



E 



*l^/ W .)(E-' b[(an 

I n=0 



2cos[(2n + 1)0„] 



7T 



2 1 ^-r. 2// -i I 7r(2n + l) 2 ^ 

This field is illustrated in figure [2j The Fourier series represents the summation of 
top-hat and saw-tooth functions. Terminating the expansion at n — n max produces 
the high frequency ripple seen in the figure. 
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3.2. Equations of motion 

The trajectory is treated classically: the electron is taken to be a point particle 
moving in prescribed external electromagnetic fields associated with the laser beam, 
and radiating energy continuously rather than in discrete jumps. When the electron 
becomes relativistic, radiation reaction plays an important role in the dynamics, and 
we incorporate it using the Landau-Lifshitz prescription [21]. To lowest order in I/7, 
where 7 is the Lorentz factor, this leads to a force / r ad that is anti-parallel to the particle 
momentum p. Writing this in terms of the unit vector in the direction of the momentum, 
p = pn one has: 

3 4 

3m 4 c 5 



/ rad = -^tIVIV (23) 



2 m2c3 2 ,OA^ 

= ~ -at—j-rfu (24) 

It can be shown [22] that the Landau-Lifshitz prescription for the classical radiation 
reaction force (the right-hand side of fl24j) does not contain h, since af oc 1/h and 
7] oc H) is the lowest-order approximation in the small parameter rjat to the physically 
acceptable solution of the exact Lorentz-Abraham-Dirac equation. As such it is valid 
up to field strengths/energies much higher than those of interest here [231 Hj. The 
force corresponds exactly to the rate of transfer of momentum from the electron to 
the photon field according to the classical formula (Q. However, the concept of a 
continuous classical trajectory fails at 77 ~ 1, where a single synchrotron photon takes 
off a significant fraction of the electron energy. Furthermore, quantum effects reduce the 
average power radiated, modifying it by the function g(rj) defined in (JT51) and plotted in 
figure lA"2l In this paper, we neglect the quantum fluctuations in the electron orbit, but 
take account of the reduction in radiated power by multiplying the radiation reaction 
force given in ( 1241) by g{rj) . The equations of motion for a particle of charge q, mass m 
are then: 

P% = -*-[E ± + frAB] 

at ^mc 

ld7 = (_q_\ g _ggWg(g) (25) 
7 dt y^mc J h 37 

where c/3 = c(7 2 — l) 1 ^ 2 / / y is the particle three speed, Ej_ is the component of E 
perpendicular to /x, and the radiation reaction term is given to lowest order in I/7. Note 
that the parameter r\ is determined by the component of the Lorentz force perpendicular 
to fx: 

7 



V 



E, 



crit 



[E ± + A Bf + U> ■ Ef /7 »] " 2 « (26) 

Given E and B as functions of position and time, equations (125]) are integrated forward 
in time using a standard fourth-order Runge-Kutta algorithm. 
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Trajectory at B=0 node 




t (laser phase) 

Figure 3. The trajectory of a particle initialized at a magnetic node in very long 
counter-propagating, circularly polarized pulses. Also plotted is the dependence of the 
parameter rj on time (measured in units of the laser phase) . 

3.3. Pair creation 

In addition to the electron's trajectory, we are interested also in the number and 
frequency of photons it radiates and, in particular, in the number of pairs created both 
by these photons as they propagate out through the laser beams and by direct trident 
pair production, according to (fT7|) . 

The number of directly created, trident pairs is easily found by adding equation 
( TT7|) to the set (I25I) . This describes the monotonic growth in time of the number of pairs 
created by intermediate virtual photons iV ± (t)| virtual . 

Pairs created by real photons are more difficult to compute. Suppose an electron 
emits a photon at time to, position xo,yo,zo with frequency v (all measured in the lab. 
frame). The total optical depth to absorption for this photon is 

mc mc /"tcsc 

T(u,to,x ,yo,z ) = at— — — / dtx(t)T±[x(t)} (27) 

n riu Jt 

where the integration is along the photon's ray path from emission to escape from the 
system at time t esc , and T± is defined in ([20]) . In general, x{t) depends on the photon's 
frequency v and its direction, which are constant along the ray path, as well as on the 
local electromagnetic field, which varies along this path. In analogy with (|26l) . % may 
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Off-node trajectory 
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Figure 4. The trajectory of a particle initialized at z = 0.5, between the nodes of 
the standing wave, in very long counter-propagating, circularly polarized pulses. The 
trajectory moves to a node of the electric field, located at z = ir/2 (indicated by the 
dotted line). Also plotted is the dependence of the parameter r] on time (measured in 
units of the laser phase) . 



also be defined as 



x(t) 



hv 



Ej_ + k A B 



(2f 



2mc 2 E ciit 

where now E± refers to the component of the electric field perpendicular to the 
propagation direction k of the photon. Therefore, assuming the photon is emitted 
in the direction of motion of the electron, from (J26l) and (l28l) we can write the photon 
frequency in terms of its x _va l ue at birth and the parameters 7 and 77 of the parent 
electron at that instant: 

w x(t )i (29) 

2mc 2 rj 

The value of \ a ^ any point on the ray path is then easily found from the local value of 
the electromagnetic fields, since both v and k are constant along this path. The total 
pair production probability per electron (via real photons) is then given by an integral 
of the pair production rate over the frequency of the emitted photons. The number of 
these pairs can also be described by a monotonically increasing function iV±(i)| rea i of 
the time t at which the pair-creating photons are emitted. This can be evaluated by 



Pair production in counter-propagating laser beams 



13 



N 



2 - 



1 - 











Off-node trajectories 



z =0.1 
z =0.5 
z =1.0 
z =2.0 




10 

t (laser phase) 



15 



20 



Figure 5. Trajectories of particles initialized at various positions in the standing 
wave in very long, counter-propagating, circularly polarized pulses. The B = node 
at z — appears to be unstable, whereas particles are attracted to the E = node at 
z = it/2. 



adding another equation to the set (|25|) : 
dN± 



dt 



v/% diV 

dx^[l-exp(-r)] (30) 



real 

Because the conversion probability, 1 — exp(— r), depends on the photon's path in 
the strong electromagnetic fields of the laser, a precise treatment requires a model of 
the spatial structure of the laser beams transverse to their direction of propagation. In 
this paper, we assume the pulses have cylindrical geometry. Taking the z-axis as the 
direction of propagation, the fields are assumed to be independent of x and y provided 
x 2 + y 2 < R 2 , where R is the radius of the cylinder, which we here set equal to one 
laser wavelength: R = 2tt, and to vanish outside this cylinder. In addition, we formally 
impose an upper limit on the length of the cylinder \z\ < -2 max , although this is not 
relevant for the trajectories discussed in this paper. The cylindrical geometry of the 
pulses enters into the computations in two ways: particle trajectories are terminated 
when they leave the cylinder, and pair creation by photons is switched off when they 
leave the cylinder. We ignore diffraction effects by treating the wave as planar within 
the cylinder, but this will have little effect on the overall results. The main motivation 
for imposing the cylindrical boundary is to take approximate quantitative account of the 
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Pair production on a circular trajectory 




23 23.5 24 
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Figure 6. The number of pairs produced during a laser period by an electron following 
a circular trajectory at the magnetic node of the standing wave, as a function of the 
intensity of (one of) the counter-propagating, circularly polarized beams. Both the 
virtual and real photon channels are shown. In addition, the parameter r\ defined in 
(|2|), which is constant on this trajectory is shown. For comparison, the number of 
pairs predicted in is shown as a dashed line. 

escape of electrons and photons from the strong field regions in the intersecting beams. 

This procedure should yield a realistic estimate of the number of pairs created 
by a single electron and the photons it radiates. However, in section 14.11 for the 
purpose of comparison we adopt an approximate method of computing the optical 
depth to absorption r that is closely similar to the method used by Bell & Kirk [I]. 
The approximation consists of replacing the integral over the photon trajectory in ( 1271) 
according to 

r<tt X T±{ X )^ — <X>7±«x» (31) 

Jt C 

The two parameters (%) and L± can be interpreted as the effective value of x averaged 
along the escaping ray path and the effective length of this path (assumed to be the 
same for each photon independent of its frequency or the position at which it is created). 

As it propagates through the laser beams, the photon encounters electromagnetic 
fields that vary in both magnitude and direction. If, for example, the photon is emitted 
almost parallel to the local electric field, the initial value of \ is much lower than the 
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Figure 7. The number of pairs produced at t = 20 (by which time this quantity has 
saturated) by an electron initialized at t — between the nodes of the standing wave, 
as a function of the intensity of (one of) the counter-propagating, circularly polarized 
beams. Only the real photon channel is shown. Initial conditions of the trajectories 
are as in figure [5] 



average value it experiences on its escape path. To allow for this, we choose the average 
value (x) so that it corresponds to a photon propagating perpendicular to an electric 
field whose strength equals the amplitude of the oscillating field at the point of emission. 
Formally, this involves summing the envelope functions (/+ and /_) introduced in ( J2TT) : 

(x) (v,t 0l x ,y ,z ) = 
hv 

2mC 2 ^ [^+(^ laserZ ° ~ ^laser^o) + /-(^lascr^O + ^laser^o)] 

so that, using ([27]) and ([ST]) , the approximate optical depth becomes 
(r) (u,to,x ,yo,zo) = 

-Q.{aT± ((x)) [j+t^laserZo ~ Wi aser to) + /-(^lascr^O + ^lascr^o)] 



(32) 



L± 



(33) 



This quantity then replaces r in ([30]) . In [1], the effective path length was estimated to 
be one laser wavelength, so that we set L± = 27rc/o;i aser . 
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4. Results 

Our results are divided into two parts. In the first, we consider counter-propagating, 
circularly polarized pulses of infinite duration. This enables us to compare with the 
results of [4]. Pair production in this case is treated using the approximation given in 
( I3T1) . In the second part, section fl~2l we consider a more realistic model, with pulses 
of finite duration. Pair production in this section is computed by integrating along the 
escape path of each photon at each timestep, and the laser fields are assumed to be 
confined to a cylindrical region. 

4-1. Circularly polarized, long duration pulses 

Bell & Kirk [4J analyzed a particularly simple trajectory: assuming two long pulses 
L — > oo, they found that at the magnetic node (B = 0), the particle settles into a 
circular orbit. This is illustrated in figure [31 The electron starts at t — at the origin, 
which is a node of the B field, with velocity directed along the positive y axis, and 
Lorentz factor 7 = 10. The intensity of a single laser pulse in this and the following 
examples is 6 x 10 23 Wcm -2 . Within less than one hundredth of a laser period the 
particle turns to move in the negative x-direction, and, after a transient lasting less 
than one half of a laser period, the trajectory relaxes to a circle in the x-y plane, with 
constant Lorentz factor 7 = 760. 

Initializing an electron at a magnetic node is convenient for analytical calculations. 
However, it selects an electron trajectory that is particularly efficient at pair production. 
An electron initialized at rest at a node of the electric field, for example, remains at rest 
indefinitely. If a starting point in between the electric and magnetic nodes is chosen, the 
electron tends to drift towards the electric node, and lose energy in the process. This is 
illustrated in figure HJ which shows the trajectory of an electron initialized at z = 0.5 (in 
units of Ai aser /27r), with the same velocity as in figure [3j As well as a drift perpendicular 
to the laser beams towards positive x and negative y, this trajectory starts to oscillate 
around the E — node at z = tc/2, and slowly converges onto it. After a sharp rise 
associated with the rapid initial acceleration experienced by the electron, the parameter 
77, that controls pair production, decreases steadily with time. 

Convergence on the E = node is a general property of the trajectories, as shown 
in figure [5] These examples demonstrate that the majority of the trajectories move 
swiftly to the vicinity of the E = mode, arriving in a zone around it of dimension less 
than one tenth of a wavelength within a single laser period. The subsequent motion is 
complex. As the electrons lose energy, their natural time and length scales decrease. 
After a long period close to such a node, this behaviour can place strong demands on 
the numerical integration method. 

The number of pairs produced by a single electron following the trajectory 
illustrated in figure [3] for a range of laser intensities is shown in figure EJ The initial 
conditions were chosen to place the particle on the circular path analyzed in [3], thus 
avoiding short-lived transients. On such a trajectory, the parameter 77 is constant. 
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t (laser phase) 

Figure 8. Sample trajectory in circularly polarized pulses of finite length, propagating 
along the z-axis. The upper panel shows position, the lower panel the Lorentz factor, 
the number iV of pairs produced per primary electron and the parameter 77 defined in 
(HI). When the pulses overlap the electric field vanishes at z = (2n + 1)tt/2, shown by 
dashed lines for n = — 2, — 1 and 0. 

The number of pairs increases linearly with time, but, for comparison with [4J, we 
select the value after one laser-period (i = to t — 2n). This figure was computed 
using the approximate description of pair production described in section 13.31 with 
L± = 2nc/u\ ascT . It is equivalent to the method used by Bell & Kirk [4] in their figure (1), 
in which, in addition, a monochromatic approximation to the synchrotron emission was 
adopted. As expected, allowing for the fact that even low energy electrons occasionally 
emit photons capable of pair producing leads to an increase in the number of pairs 
produced at low laser intensities. As a result, the pairs produced via virtual photons 
are overwhelmed by those produced by real photons, except at / < 10 23 Wcm -2 , where 
the production rate is in any case very low. An additional factor contributing to the 
higher rate of pair production found in these calculations when compared to jl] is the 
inclusion of quantum corrections to the radiation reaction force via the factor g(r)) in 
( |25|) . This leads to larger values of both 7 and r\ for given laser intensity. 

If the electron is initialized in between the nodes, pair production is reduced. This 
is shown in figure [7J In contrast with the circular trajectory at the B = node, the 
number of pairs produced when the particle is initialized off-node does not increase 
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t (laser phase) 

Figure 9. Sample trajectory in finite-length pulses with aligned, linear polarization. 

linearly with time, but saturates after a few laser periods. Consequently, figure [7] shows 
the time- asymptotic number of pairs, for the four initial conditions shown in figure 

The curve for zq = 2.0 demonstrates that, for circular polarization, electrons 
initialized close to the E = node move towards it rapidly without emitting as 
many pairs as those electrons initialized close to the B = node. Electrons with 
an unfavourable initial position do not contribute substantially to pair production. 
However, this negative effect is offset by the positive modifications illustrated in figure El 
The general conclusion of [4J that an ensemble of N electrons produces approximately 
iV secondary pairs when the intensity of circularly polarized counter-propagating laser 
beams into which they are injected reaches 10 24 Wcm -2 remains unchanged. 

4-2. Finite duration pulses in cylindrical geometry 

The results of the previous section show that the position of a particle within a circularly 
polarized, long pulse train is important in determining the number of pairs it produces. 
In reality, however, pulse lengths are limited to a few laser wavelengths. As the pulses 
approach each other they pick up electrons and carry them forward on the front edge 
of the pulse until interaction starts. During interaction, the electron trajectory remains 
for an extended period at approximately constant z. 

In the case of circular polarization, the particle stays close to a node of the electric 
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field, but in the case of linear polarization, it samples a much larger region. In general 
terms, this is due to the characteristic figure-of-eight oscillation of an electron in a 
linearly polarized wave during which the electron oscillates in the direction parallel to 
beam propagation as well as in the perpendicular direction. In contrast, in a circularly 
polarized wave the electron can move smoothly, without parallel oscillation, towards 
the attracting E = node. This behaviour is illustrated in figures [8] and [91 in which 
the laser intensities and laser pulse shapes are the same but the polarizations differ. 
In each case, pulses of peak intensity J24 = 0.6 with length L = IOit and thickness of 
the pulse edge A = 107r/3 were chosen. The trajectory was initialized at t = —100, 
when the pulses are well separated, at the position x — y — 0, z — —50, moving in 
the positive z direction with Lorentz factor 10. The particle coasts at constant 7 until 
t ~ —30, when it encounters the leading edge of the leftward propagating wave. This 
wave picks it up, reversing its motion along z and accelerating it until the pulses begin 
to interact at £ ~ 0. In the circularly polarized case, the particle immediately drops 
onto the E = node at z = — n/2 (shown as a dashed line), the Lorentz factor falls and 
few pairs are created (~ 10~ 7 per primary particle). However, in the linearly polarized 
case, the particle performs substantial excursions in z, settling only briefly onto the 
E = node at z = — 37r/2. During this time, the Lorentz factor fluctuates rapidly, with 
only a slight decrease in its average value. As a result, a relatively large number (0.3 
per primary particle) of pairs are created. The trajectories terminate when the particle 
exits the cylinder containing the laser pulses at roughly t = 50 for circular polarization 
and t = 65 for linear polarization. 

During the period when the pulses overlap, the electron trajectory is very complex, 
especially in the case of aligned linear polarisation. This makes it quite sensitive to 
the accuracy demanded of the integration algorithm, in the sense that a more stringent 
accuracy requirement leads to a trajectory that remains converged to larger times. For 
the trajectory shown in figure [9] we find that an error per step of less than 5 = 10~ 4 
in position and in the angles (in radians) used to describe the unit vector /x yields a 
converged trajectory up to t = 20, whereas 5 = 10~ 6 extends this to t = 35. (The 
example presented in figure [9] was computed with 5 = 10~ 8 .) However, in both cases 
pair creation ceases well before accuracy is lost. We conclude that counter-propagating, 
circularly polarized beams are intrinsically inefficient at pair creation, at least when the 
fields rotate in the same direction, and do not consider them further. 

Trajectories in linearly polarized pulses of finite length, contained in the cylindrical 
volume described above are sensitive to the point at which they are first picked up by a 
pulse, and to the time they subsequently spend at its leading edge, before they encounter 
the counter-propagating pulse. To take account of this, we present, in figures [101 EH 
and [121 results obtained in pulses with L = 10n and A = 107r/3 by varying the initial 
x, y and z positions of the particle, whilst keeping the other initial conditions fixed. 
In each case, 3 x 10 4 trajectories were computed, using randomly selected values of x 
and y uniformly distributed over the cross-section of the cylinder, random values of z 
distributed uniformly over the range —50 to 0, and random values of the pulse intensity, 
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Figure 10. Pair production rates for a range of initial particle positions and for aligned 
linear polarization. Each pixel corresponds to a small range of laser intensity ("flux", 
shown on the x-axis) and of the number -ZV rca i of pairs produced via the channel that 
involves an intermediate real photon (shown on the y-axis). The colour coding of a 
pixel indicates the logarithm (to base 10) of the probability that, given a flux in this 
interval, an electron with randomly chosen initial conditions (see text) will produce the 
corresponding number A^ rea i of secondary pairs. Roughly 80% of trajectories produce 
no pairs and so are not shown. At a laser flux of 10 24 Wcm~ , about 10% of the 
particles reach the interaction region and the majority of these produce on average 
roughly one secondary pair. 

distributed uniformly in logarithm over the range 23 < log 10 (^24) < 24. The initial 
Lorentz factor is set to 10 a001 , corresponding to an electron energy of about IkeV, and 
the initial velocity is in the positive z-direction. In each of the three cases, roughly 80 % 
of the trajectories exit the cylinder before reaching the region where the laser beams 
interact, and so fail to produce pairs. The remaining trajectories were accumulated in 
60 bins in intensity and 60 bins in pair yield, and are displayed as colour contour plots. 

The aim in figures [TQl [11] and [12] is to test the sensitivity of pair production to 
polarization and pulse shape. In figure HH "crossed polarization" , the electric field of 
one beam is in the x-direction whereas the field of the other beam is in the ^-direction. 
In figures [TQl and [121 the electric field of each pulse is in the x-direction. Figure [TQl and ITT1 
present results obtained using sinusoidal waves (modulated by the envelope function), 
whereas in figure [12] high harmonics are included in one of the waves, as illustrated in 
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Figure 11. Pair production rates for a range of initial particle positions, as in figure[T(71 
but for crossed linear polarization 

figure 

The possible advantage of the crossed configuration is that only at isolated points 
and times is the total electric field zero. When considering circular polarization in section 
14.11 it was found that electrons congregated at the E = node and pair production was 
dramatically reduced. By eliminating the E = node this possibility is removed, which 
suggests that pair production might be increased. However, comparison of figures [10] 
and [TT] shows that "aligned polarization" with both electric fields oscillating in the 
same direction performs nearly as well as crossed polarization - possibly because the 
maximum electric field is larger by a/2. In all three cases, (figures dUJ dU and [T2l) those 
trajectories that produce fewer than 10 -10 pairs per primary electron are omitted. The 
remaining trajectories, show a large range in the number of pairs created, particuarly 
for aligned polarization at low laser intensity. However, most of the electrons sit on 
a band stretching from iV rea i rs 10~ 8 at an intensity of 10 23 Wcm~ 2 to iV rea i « 1 at 
10 24 Wcrn -2 . This implies that pair production is effective for those particles that reach 
the interaction region, although only a minority of electrons do so if, as we assume here, 
they are injected uniformly over the beam cross-section just before the pulses meet. 

In figure [T2J the hypothesis is tested that pair production can be enhanced if the 
waveform contains high-frequency harmonics. At first sight, this is not implausible, 
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Linear polarization, reflected wave 
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Figure 12. Pair production rates for a range of initial particle positions as in figure ITOT 
but for aligned linear polarization with the leftward propagating wave modified to 
represent a reflected wave (with n ma x = 7, see (f22]l ). 

because of the following argument: As discussed in |4], radiation losses prevent electrons 
reaching the Lorentz factor 7 max = eA max /mc 2 , where A max is the maximum value of 
the electromagnetic vector potential. They do this by aligning the particle velocity and 
the accelerating force exerted by the field, and, therefore, not only reduce the Lorentz 
factor below ev4 max /mc 2 , but also reduce the perpendicular component of the field and 
the parameter rj that controls both pair production and radiation losses. Prolific pair 
production would be much easier to achieve if radiation losses could be avoided whilst 
the electron is being accelerated. This suggests that if one of the counter-propagating 
waves contained abrupt changes in the electric or magnetic field it might be possible 
for the other wave to accelerate the electron in the direction of laser propagation to a 
high Lorentz factor but relatively small rj. Subsequently, an abrupt change in the field 
direction caused by the first wave could lead to a large value of r] before radiation losses 
set in. 

Such a situation might arise naturally when a sinusoidal wave is incident on a 
solid target since, as shown by [T7], the reflected wave contains harmonics and abrupt 
changes in the electric field. According to PIC simulations [18) . at an intensity of around 
10 21 Wcm -2 , the reflected waveform is generically similar to that in figure [2j The 
results presented in figure [12] are computed using this waveform for one of the waves 
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while keeping the total intensity unchanged. However, at least in this example, our 
hypothesis is not confirmed: the presence of harmonics appears to make little difference 
and electrons produce roughly the same numbers of pairs as in sinusoidal waves with 
aligned linear polarization. 

We conclude from figures fTUl [TT1 and fT2l that pair production is relatively insensitive 
to the orientation of the planes of linear polarization and to the waveform. While these 
results provide no indication that we can further increase pair production by a careful 
choice of waveform, they reinforce our general conclusion that strong pair production 
takes place whenever counter-propagating or reflected beams reach intensities around 
10 24 Wcm" 2 . 

5. Conclusions 

In the above, we examine the possibility that an electron moving in counter-propagating 
laser beams of intensity ~ 10 24 Wcm -2 could produce a substantial number of electron- 
positron pairs solely as a result of interaction with the laser fields. This is done using 
an accurate treatment of the physical processes, and embedding these in a numerical 
analysis of classical trajectories in realistic field configurations. 

In the case of the photon production rate, we lift the monochromatic approximation 
used in [I] and take full account of quantum effects in the weak-field, quasi-stationary 
approximation (that is well-justified in the regime of interest). The result is shown 
in figure where we examine a circular trajectory at the B = node of counter- 
propagating, monochromatic, circularly polarized beams and compare our results with 
[3]. Because even relatively low energy electrons produce a few hard photons that 
are subsequently capable of pair production, pairs start to appear at even lower 
intensities than previously predicted. The process of direct production via virtual 
intermediate photons, which was predicted to dominate over production via real photons 
for intensities less than 3 x 10 23 Wcm -2 is swamped by this effect for laser intensities 
exceeding 10 23 Wcm~ 2 . At lower intensities, the virtual photon channel still dominates, 
but the overall rate of pair production is very small. At higher intensities, where 
the parameter t] is of the order of unity, classical synchrotron theory predicts the 
emission of photons whose energy exceeds that of the radiating electron. When this 
is corrected, (see figures [H and IA2j) . not only is the rate of emission of hard photons 
reduced, but the radiation reaction force acting on the electron is diminished. Radiation 
reaction significantly inhibits pair production on a circular trajectory, so that these two 
corrections have opposite effects on the pair production rate. Figure [6] shows that the 
overall effect is a slight enhancement of the pair production rate — we find the threshold 
above which an electron on this circular trajectory produces more than one pair per laser 
period is crossed at about 5 x 10 23 Wcm -2 , compared to 10 24 Wcm~ 2 found in [3]. 

To take account of more realistic trajectories, we consider pulses of finite length 
and finite cross-section, with both linear and circular polarization, and inject electrons 
at a complete range of phases and positions. In addition to sinusoidal waveforms, we 
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examine a laser pulse with harmonics, such as might be generated by reflection by an 
overdense plasma. The electron-laser interaction at the B = node considered in [1] 
is a special case, not least because the trajectory is unstable and the electrons migrate 
towards the E = node if the laser beams are circularly polarized. Nevertheless, 
under a wide range of laser pulse-shapes and polarizations with injection at random 
phases, we find that the conclusion is essentially unchanged: as the intensity of counter- 
propagating beams approaches 10 24 Wcm -2 , the number of secondary pairs produced 
becomes roughly equal to the number of primary electrons interacting with the beams. 
Since each electron and positron so produced interacts further with the laser beams to 
produce further pairs, an avalanche of pair production is possible. 

Further work is needed, particularly to remove the assumption that electron 
radiation losses are continuous, and to follow the interaction of primary and secondary 
particles with the complicated field structures produced by laser-solid interactions. 
However, our results strongly suggest that it may be possible to produce substantial 
pair plasmas with the next generation of extremely high-power lasers. 
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Appendix A. Synchrotron emissivity 

The function F(j],x) is given in equations (2.5a. . .f) in [19] : 

F(v, x) = E mxh) My) (A.i) 

where 

y = [*n(v-2x)} (A ' 2) 

is the quantum equivalent of the classical parameter (photon frequency /characteristic 
synchrotron frequency). The Fi are all positive, provided that the electron energy 
exceeds that of the photon, r\ > 2%: 

F 1 (x) = i + (i- x y 2 

F 2 (x) = 2{l-x)- 1 

F 3 (x) = x 2 (1 - xY 2 (A.3) 



The Jj are also positive: 
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Figure Al. The synchrotron functions Ji defined in (|A.4|) 
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For y <C 1 they are well approximated by 
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These functions are plotted in figure | 

The function g(rj) defined in (|T5|) that governs the total radiated energy can be 
written 



(A.6) 
(A.7) 



g(v) = ^ / 



2y 2 K 5/3 (y) 36 V 2 y 3 K 2/3 (y) 



87T JO 

This function is plotted in figure [A2 
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Figure A2. The function 17(77) defined in (fT5)) that describes the quantum corrections 
to the total power in synchrotron radiation 



A closely related alternative approximation to the synchrotron function is given by 



1251: 



4x ; 



-yK 2/3 (y) + (l-^ y J™ dtK 5/3 (t) 



(A.9) 



where y is defined in (1A.2D . 



Appendix B. Pair-production via virtual photons 

The function £l(rj) can be evaluated starting from the expression given in [T9 
(equation A39): 



(■OO 

n"(x) =4vr / &uK 2 (Axu)K 2 ,z(l/u) 
Jo 



(B.l) 



Mathematica gives an explicit expression for this integral in terms of the 
Meijer G function, and integrates it twice (using the boundary conditions Q'(oo) = 
fi(oo) = 0) to obtain 
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